Dissociation dynamics of resonantly coupled bose-fermi mixtures in an optical lattice 
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We consider the photodissociation of ground-state bosonic molecules trapped in an optical lattice 
potential into a two-component gas of fermionic atoms. The system is assumed to be described by 
a single-band resonantly-coupled Bose-Fermi Hubbard model. We show that in the strong fermion- 
fermion interaction limit the dissociation dynamics is governed by a spin-boson lattice Hamiltonian. 
In the framework of a mean-field analysis based on a generalized Gutzwiller ansatz, we then examine 
the crossover of the dissociation from a regime of independent single-site dynamics to a regime of 
cooperative dynamics as the molecular tunneling increases. We also show that in the limits of weak 
and strong intersite tunneling the mean-field solutions agree well with the results from the quantum 
optical Jaynes-Cummings and Tavis-Cummings models, respectively. Finally, we identify two types 
of self-trapping transitions, a coherent and an incoherent one, depending on the ratio of the repulsive 
molecule-molecule interaction strength to molecular tunneling. 

PACS numbers: 03.75.Ss, 05.30.Fk, 32.80Pj, 67.60.-g 



I. INTRODUCTION 

The formation of molecules by Feshbach resonance m 
and by two-photon Raman photoassociation from ul- 
tracold atomic gases trapped in an optical lattice po- 
tential offers remarkable opportunities to study some of 
the long-standing questions of condensed matter physics 
in systems under exquisite control, and is also of much 
interest for quantum information science 0, 0| . Exter- 
nal control parameters such as laser and magnetic field 
strengths can be controlled precisely in real time, al- 
lowing the experimental investigation of the dynamics 
of strongly correlated systems from the adiabatic to the 
sudden regime, the study of the formation of long-range 
order and defects as the system transits across quantum 
critical points, and the engineering of exotic many-body 
states starting from various quantum phases ^ ^] . 

The creation of molecules from cold atoms in optical 
lattices has been demonstrated experimentally by sev- 
eral groups. Applying two-photon Raman lasers on *^Rb 
atoms in a 3-dimensional (3D) lattice, T. Rom et al. 
succeeded in controlling both the internal and center- 
of-mass states of the molecules by changing the detuning 
of the photoassociation laser fields 0- Unfortunately, 
the molecular decay arising from the spontaneous emis- 
sion induced by Raman laser photons exceeded the pho- 
toassociation coupling, so that this experiment could not 
observe coherent Rabi oscillations between atoms and 
molecules. The first lattice experiments exploiting the 
passage of fermions ('^^K atoms) across a Feshbach res- 
onance were performed by M. Kohl et al. 0- They ob- 
served the occupation of atoms in higher Bloch bands 
near the resonance, indicating that the inclusion of inter- 
band coupling is crucial for the description of a broad 
Feshbach resonance in optical lattices H . They also cre- 
ated Feshbach molecules from fermionic atoms in a three- 
dimensional optical lattice ^lO.J . and found that the mea- 
sured binding energy of the molecules is in good agree- 
ment with the theoretical predictions of single-channel 



models of atoms in a tight trap 0, . We also men- 
tion a recent experiment by G. Tialhammer et al. |l3j| . 
who created long-lived *^Rb2 molecules, with lifetimes of 
up to 0.7s, in a three-dimensional optical lattice via a 
narrow Feshbach resonance. 

In a recent experiment, C. Ryu et al. generated a 
quantum degenerate gas of molecules via photoassocia- 
tion of ^^Rb atoms, and observed a clear coherent atom- 
molecule oscillation in an optical lattice, with Rabi fre- 
quency « 4.13kIIz 14]. In this experiment the atoms 
were initially in the Mott-insulator (MI) phase. For the 
duration of the photoassociation pulse each atom was 
trapped in the ground state of a single lattice site with 
trapping frequency ^lattice ~ 27r x 27kIIz, and the atom- 
atom interaction strength and atomic intersite tunneling 
were of the order of ~ lOkHz and ^ lOHz, respectively. 
In general, the initial atomic state in that experiment 
can be in either the MI or the superfluid (SF) phase, de- 
pending on the depth of the optical lattice. Combined 
with the interplay between two-body interactions, inter- 
site tunneling, and atom-molecule coupling, this results 
in the photoassociation process exhibiting a very rich dy- 
namics. 

Motivated by these experimental advances, this pa- 
per considers a gas of two-component fermionic atoms 
coupled to bosonic molecules via photoassociation in an 
optical lattice. The system consists initially of bosonic 
molecules only, assumed to be in a ground state corre- 
sponding to either a MI phase or a SF phase, and a cw 
photoassociation beam is switched on instantaneously at 
t — Q. Our goal is to study the subsequent dissociation 
dynamics of the bosonic molecules into fermionic atoms. 
We examine the crossover of the dynamics from an inde- 
pendent single-site regime to a cooperative regime. We 
also show th e app earance of two types of self-trapping 
transitions |l5l llq . which we characterize as coher- 
ent and incoherent, as the ratio of repulsive molecule- 
molecule interaction strength to molecular tunneling is 
varied. 
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Section II discusses our model and formulates it as 
a single-band resonantly coupled Bose-Fermi Hubbard 
model. We show that in this model the dissociation dy- 
namics can be described by a spin-boson lattice Hamil- 
tonian in terms of pseudo-spin formalism. In this section 
we also give a formulation of mean field approximation 
based on a generalized Gutzwiller (GW) ansatz. Section 
III discusses the dissociation dynamics in several limit- 
ing cases. Specifically, we show that in the weak and 
strong molecular tunneling limits, respectively, it can be 
described in terms of generalized versions of the Jaynes- 
Cummings and Tavis-Cummings models of quantum op- 
tics. Section IV then presents numerical results of the 
GW mean- field dynamics, and section VI is a summary 
and outlook. 

II. MODEL 
A. Spin-Boson lattice Hamiltonian 

We consider at zero temperature a gas of fermionic 
atoms of mass m/ and spin c =t, i trapped in an optical 
lattice potential. The fermionic atoms can be coherently 
combined into bosonic molecules of mass mf, = 2m/ via 
two-photon Raman photoassociation. The lattice lasers 
are adjusted in such a way that the tight binding approx- 
imation is well justified. Since in lattice photoassocia- 
tion experiments the atom-molecule interaction energy is 
typically small compared to the trapping energies of the 
atoms and molecules 0, 0| (note that this is not so for 
the broad Feshbach resonance of atoms HITqI) the 
system can be described by the sin gle- band resonantly 
coupled Bose-Fermi Hubbard model [TtL lll [T9ll20l| 

HBF = Hf + Hb + Vbf, (1) 

where 

Hf - -hJf CfUoo + B.C.) + hcvj J2 

{ij)a i<y 

i 

Hb = -hJbY{brbj+}i.c.) + h{uJd + uJb)Yn'l 

(ij) * 
i 

Vbf = hUbfY.nln': + hgJ2ifyiik + ii.c.). (4) 

ia i 

Here fi^ and 6; are annihilation operators for the 
fermionic atoms of spin a and bosonic molecules at the 
i-th site, respectively. The corresponding number op- 
erators n{„ and have eigenvalues n{„ and n^. The 
conserved quantity 

Nr^'^Y.^+Y.^L (5) 

i ia 



is the total number of fermions in the system. 

The Hamiltonians Hf and Hb describe standard Hub- 
bard models for fermions and bosons, respectively, and 
Vbf describes the interactions between fermions and 
bosons. The terms proportional to J^, where a = f,b, 
account for the tunneling of particles between nearest 
neighbor sites denoted by (ij). The single-particle center- 
of-mass energies are htUa, and the effective detuning be- 
tween the Raman lasers and the difference in internal 
energies of the atoms and molecules is hcod- The on- 
site fermion-fermion, boson-boson, and boson-fermion 
collisions are described by the interaction strengths Uf, 
Ub{> 0), and Ubf, respectively. Finally, the term propor- 
tional to g in Vbf describes the conversion of fermionic 
atoms into bosonic molecules, and vice versa, via pho- 
toassociation. 

We concentrate on the case where the system con- 
sists initially of bosonic molecules only, and the corre- 
sponding initial state is the ground state of the Bose- 
Hubbard Hamiltonian Hb- It is well-known that the 
ground state of this model describes either a MI phase 
or a SF phase, depending on the ratio Ub/zJb, where z 
is a number of nearest neighbor sites, for a fixed number 

of bosons m mm mil El. 

In the following we consider a strongly confining regime 
\Uf \ ^ zJf. In this regime, Fermi pair states, that is, 
pairs of fermions with spins up and down occupying the 
same site, see upper diagrams in Fig. ^ have an energy 
separation Uf compared to unpaired states such as in the 
lower diagram. This allows us to treat the fermionic tun- 
neling that couples the paired and unpaired states per- 
turbatively. Since in addition the photodissociation of a 
molecule always creates a Fermi pair at a given lattice 
site, the unpaired states can then be adiabatically elim- 
inated and integrated out. This results in the tunneling 
of Fermi pairs via those virtual states (Fig. QJ, and the 
Hamiltonian Hf is mapped on a pseudo-spin- 1/2 system 
(XXZ model) by second-order degenerate perturbation 
theory [11 113, 

* (ij) 

(6) 

Here cug = ojf + Uf/2 and Js = AJ'j/Uf, and we have 
introduced the pseudo-spin operators describing the cre- 
ation and annihilation of a Fermi pair 

and 

s/. = ^(4+4-l) (8) 

at the i-th site ^28.] . This mapping amounts to describ- 
ing fermionic pairs as "effective two-level atoms" whose 
upper level corresponds to a pair of fermions, and lower 
level the absence of such a pair. Thus the presence of 
paired states at the z-th site is described by the spin- 
up state \sj^ = 1/2) and its absence by the spin-down 
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FIG. 1: (Color online) Schematics of the tunneling of a Fermi 
pair in the limit Uf S> zjf. Filled and open circles present 
fermions with different spins. 

state |s^j = —1/2). The pseudo-spin operators obey the 
commutation relations 

In terms of these operators the interaction Hamiltonian 
Vbf can be re-expressed as 

Vsb = nutf ^(2J^, + i)h1 + ngY^{b\s-^ + s+S,). (9) 

i i 

We thus conclude that in the strongly confining regime 
the dissociation of initially bosonic molecular states into 
fermionic atoms is governed by the spin-boson lattice 
Hamiltonian 

HsB^Hs + Hb + Vsb. (10) 

This effective Hamiltonian is valid for both of attractive 
and repulsive fermion-fcrmion interaction cases for the 
problem at hand. 

B. Gutzwiller mean-field ansatz 




FIG. 2: Superfluid density |$|'^ of the molecular ground state 
as a function of Ub/ zJb. The MI-SF critical point is at about 
5.83. 



The time-dependent Gutzwiller approach has previously 
been used in the studies of several ultracold systems on 
lattices, such as bosonic atomic-molecular gases 30], 
disorder jsj , inhomogeneous |3j| Bose-Hubbard systems 
and non-resonant Bose- Fermi mixtures (33j |. In all cases 
the time dependence of the system parameters was as- 
sumed to be slow enough that the system is allowed to 
evolve adiabatically. 

In the present case, the ground state of the bosonic 
molecules in the absence of photoassociation is given by 
the stationary solution of Ea. H12|l with f^\/2 ~ ^ ^^'^ 
g = Q. One of the merits of the Gutzwiller ansatz is 
its prediction of the MI-SF phase transition for a com- 
mensurate filling factor. For a fixed number of bosonic 
molecules, the ground state is determined by the single 
parameter Uh/zJ^ and the corresponding configuration 

of iff' _i/2\ is uniform over all lattice sites. At Uh = 0, 
corresponding to the deep SF regime, the ground state 
is a coherent state and the variation parameters for the 
system for the filling factor v are given by 



In order to proceed we now introduce a Gutzwiller- 
type (GW) variational ansatz [22, 113 for the many-body 
wave function as 



= n 



oo 



fl:lsAi)\n,)b®\s^) 



fp 



1/2 

^ E 

1=1 \ni=0si=-l/2 

(11) 

where Ng is the number of lattice sites, and |ni){, and 
\si) fp refer to number states of the bosons and Fermi pair 
at the i-th site, respectively. The normalization condi- 
tion J2ni Si = 1 is satisfied at every single site. 
From the time-dependent variational principle, the time 
evolution of the variational parameters fn^si is obtained 
via the minimization 



d 



d 



*W = 0. (12) 



Ai) _ -./2 (V^)"' 

■'ni,-l/2 — ^ 



(13) 



leading to Poissonian number statistics at each site. For a 
finite Ub the ground-state number statistics of molecules 
becomes subpoissonian, and in the MI phase, Ub/zJi, > 

5.83 for V — becomes a Fock state f!^^ _i/2 — ^ni,v 
Fig. 121 shows the superfluid density per site |$p of the 
molecular ground state, as a function of Ub/ zJjj for the 
commensurate case v = 1, where $ = {bi). 

It can be shown from the minimization of Eq. H12I) that 

the resulting variational parameters {fn},si} remain uni- 
form over the entire lattice when starting from a uni- 
form configuration. This is the case for the dissocia- 
tion of ground-state bosonic molecules. In that case the 
equations of motion for the site-independent parameters 
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{fn.s} are found to be 

i = hn^sfn,s + g [Vnfn~^l,s+1 + ^ n + lfn+l,s-l] 

+ ^[A*/„,.+i+/„,.-iA-2sAf (14) 

where ~ 6n + {Ub/2)n{n — 1) + Ubfn{2s + 1) and 
S = LOd + ujb — 2(jjs. Here, we have used the average 
number conservation 

= 2A^s|e(" + '^ + ^)|/»,^I'| (15) 

and Nra is the initial number of molecules. We have 
also introduced the superfluid order parameters of the 
molecular bosons, $, and the atomic Fermi pairs. A, and 
the magnetization M of the pseudo-spin as 

n.s 

A(<) EE {Mk) = {Sj,)=Y.f:jn,s+l, 

71. S 

M{t) EE 

n,s 

Note that in Eq. H14|l . the tunneling of both bosons and 
Fermi pairs occurs only via their own order parameters, 
and fluctuations about their mean values are ignored. In 
this sense, the GW ansatz can be thought of as a kind of 
mean-field approximation. 



III. LIMITING CASES 

Before turning to the presentation of a numerical study 
of Eq. H14|l , we first discuss analytical results for the disso- 
ciation dynamics in several limiting cases that provide us 
with some useful intuitive understanding of its behavior. 



A. Generalized Jaynes-Cummings dynamics : 

zJb, zJs < g 

In the limit Jb, Js — > 0, the spin-boson Hamilto- 
nian Hsb reduces to a generalized version of Jaynes- 
Cummings (J-C) model |3J| at each lattice site. This 
situation can be solved exactly |2^. The average num- 
ber of molecules per site, n^, is 

n,(t)=5:p(n)^cos2(^), (16) 
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FIG. 3: Time evolution of the average molecule number per 
site rii, according to the Jaynes-Cummings dynamics for ini- 
tially (a) a number state corresponding to the Mott-Insulator 
phase and (b) a coherent state. Both situations are for 
(7b/ = and exact resonance 5 = 0, and for u — 1. The 
dashed line is the GW mean-field solution for zJt/g — 0.1. 

where the Rabi frequency TZn, 

Tin = ,J{S+ (Ub - 2Ubf)nY + ^g\n + 1), (17) 

depends on the number n of molecules. In Eq. Hlt)|l . 
■p(n) = |/n.-i/2(^ = 0)P is the initial number statistics of 
the molecules at each lattice site. We note that the GW 
mean-field equation (|14(l has the same exact solution and 
thus describes the correct dynamics in this limit. 

The time evolution of at exact resonance, (5 = for 
a molecular field initially in the MI phase {Uh ^ zJh) 
and the superfluid phase (Uh = 0) is shown in Fig. 13 (a) 
and Fig. O (b), respectively, for the filling factor = 
Nm/Ns = 1. In this example we have taken Ubf — 0. 
Note also that the on-site interaction between atoms and 
molecules plays no role in the dynamics starting from 
the MI state for = 1. The dashed hue in Fig. H (b), 
obtained from Eq. (|14|l for small but finite tunneling, 
zJb/g = 0.1, illustrates the small deviation from the 
Jaynes-Cummings solution Hlt)|) in that case. 

These results demonstrate that the observation of the 
number of atoms, which is the difference between the to- 
tal number of fermions and twice the molecule number 
and is easier to observe than that of number of molecules, 
can reveal the precise number statistics of the molecules 
and hence their ground-state properties. A similar argu- 
ment was previously given by K. M0lmer in the context 
of bosonic atom-molecule gases in an optical lattice [s^ . 



B. Generalized Tavis-Cummings dynamics : 

zJb > g and Ub/zJt < 1 

As shown in Fig. [3 for Ut/zJb < 1 the initial molecular 
ground state contains a large condensate fraction or high 
occupation of the zero-momentum state. Moreover, the 
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application of a weak photoassociation coupling, z ^ 
presumably preserves the coherent character of the 
molecular condensate. Owing to Bose enhancement, in 
that limit the dissociation dynamics is therefore expected 
to be dominated by the condensate of bosonic molecules. 

Assuming Js = for simplicity, the momentum repre- 
sentation of the Hamiltonian H1(J|I is 



SB 



E 



s^,+H.c. 



ki,k2 i 

1^ E E ^k^+k^M+k,J}i^^ibkJ>ki (18) 



where bk is the annihilation operator of the bosonic 
molecules with momentum k and is the coordinate of 
site i. The single-molecule frequency for a regular lattice 
of lattice constant a in d dimensions is given by 



= (5 — 2 Jf, cos(fcja). 



In the extreme superfluid case, the non-condensed 
molecule fraction is negligible and a single-mode approx- 
imation that accounts only for the zero-momentum state 
is appropriate. The replacement bk bo, together with 
the approximate number conservation 



(19) 



where no ~ ^o^o, simplifies then the Hamiltonian l|18|) 
to a generalized version of the Tavis-Cummings (T-C) 
model of quantum optics [s^ 



hg 



bisj 



H.c. 



nuo 

2iV. ' 



(20) 



Here Sj — J2i ^fi ^ collective spin operator over the 
entire lattice, and the effective detuning and effective in- 
teraction strength are now 

6o ^ 6 ~ zJb + 2iyUbf, Uo^Ub-AUbf, 

where we have neglected a term of 0{l/Ns). As a re- 
sult of the absence of the phase factor e*'"'"' for the zero- 
momentum condensate molecules in the conversion term 
of Hq, the system undergoes cooperative oscillations be- 
tween atoms and molecules. Similar models have been 
applied to the atom-molecule dynamics across the BCS- 
BEC crossover and the dissociation dynamics of 

a molecular BEG .39] in free space. 

Although the Hamiltonian H2U|I can be solved numer- 
ically by direct diagonalization, it is useful to consider 
an approximate analytical solution valid in the semiclas- 
sical limit ^| . By treating the molecule number opera- 
tor no — > no classically and taking into account number 



conservation (|19|l and energy conservation, we can recast 
the Tavis-Cummings dynamics in the form of a Newton 
equation 



rno(t) 



-V{no), 



where the potential V{n) is given by 



(21) 



V{n) 



+1 TV + 2A^/" 



1 



2.9 {l-iy)+6o + 



2 4 



(22) 



In the case Uq = 0, the potential is a cubic function 
of n, and no(t) is given in terms of the Jacobian elliptic 
function sn(0; k) [41| as 



"o(^) 

where 
n± 



1 - (1 - n+)sn2 V(l - n^)iygt- , (23) 



1 

2^ 




1^ + 1 



± 




and k — ^(1 — n+)/(l — n_). 

For Sq — 0, the classical trajectory approaches an un- 
stable equilibrium point at no = when v = I 38]. It 
can be shown from Eq. H23|l that v 1 excludes the 
presence of unstable extremum. For large detunings, 
Sa/2g ^ y/V, 1, the oscillations become sinusoidal 



"o(^) 



-'I 



sm 



where 



n 



^5l+A{l + v)g^ 



(24) 



(25) 



Figure0]shows the time evolution of n?, for 5{)/ g — —1, 
Ub = Ubf = 0, and v = \. The solid and dashed lines 
are the semiclassical and quantum-mechanical Tavis- 
Cummings solutions, respectively, for N„i — 2000 and 
nb — Uq. The GW mean-field solution (dot-dashed line 
for zJb/g = 20) is seen to agree quantitatively with 
the Tavis-Cummings dynamics in the limit zJb S> g for 
Ub/zJb < 1. 



IV. GUTZWILLER MEAN-FIELD DYNAMICS 

The preceding section illustrated in limiting cases the 
crossover of the dissociation dynamics from a regime of 
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FIG. 4: (Color online) Time evolution of for 5o/g = —1, 
Ub = Ubf = 0, and u = 1. The solid and dashed lines corre- 
spond to the semiclassical and quantum {Nm = 2000) solu- 
tions of the Tavis-Cummings dynamics, respectively, and the 
dot-dashed line is the GW mean- field solution for zJb/g — 20. 



independent single-site dynamics to a regime of coopera- 
tive dynamics as the molecular tunneling is increased. In 
this section we examine the characteristics of the disso- 
ciation in these regimes in more detail within the frame- 
work of a Gutzwiller ansatz. We also demonstrate the 
appearance of two types of self-trapping transitions, de- 
pending on the ratio Ub/zJb- We restrict our discussion 
to the case Jg — Ubf — and to the commensurate case 
J/ = 1 for simplicity. The effects of the XXZ term propor- 
tional to Js in the Hamiltonian are briefly discussed 
in the last subsection. 



A. Crossover between single-site and cooperative 
dynamics 

Figure |S1 shows a schematic phase diagram summariz- 
ing the characteristics of the GW dynamics in the Uz/g 
- zJb/g plane. For zJf, 3> g and Ub/zJb ^ 1, the dynam- 
ics is dominated by the molecular condensate and hence 
cooperative. In the opposite limit, zJb <SC g, the dynam- 
ics on a given lattice site becomes increasingly indepen- 
dent of the other sites. Note that the thick gray lines 
in the figure do not indicate sharp boundaries between 
these regimes. The dashed line marks the separation be- 
tween initial molecular ground states in the the MI and 
SF phases. Note that for a molecular field initially in 
the MI phase [Ub/zJb ^ 5.83), the full Rabi oscillation 
of Fig. is always recovered at every single site since 
{bi) = in that case. 

Figure |H| shows time evolution of the average number 
Ub of molecules (solid line), the superfluid density |<i>p 
of bosons (dashed hne), and the order parameter |Ap of 
Fermi pairs (dot-dashed hne) per site for (a) zJb/g = 0.1, 
(b) zJb/g — 1, and (c) zJb/g = 10 and for 5 — zJb and 
a fixed Ub/zJb = 0.1. The corresponding probabilities 
|/n,sP for these three regimes are shown in Fig. [7| The 



FIG. 5: (Color online) Phase diagram of the GW mean-field 
dynamics. The dashed line indicates the MI-SF phase transi- 
tion for the initial molecular ground state. 



initial molecular ground state is the same in all cases and 
corresponds to an almost coherent state. 

In the case of weak tunneling, Fig.|Hl(a), the system dy- 
namics resembles the Jaynes-Cummings limit, see Fig. |31 
(b): pairs of amplitudes, fn-1,1/2 and /„,_i/2, are cou- 
pled to each other but uncoupled to the other amplitudes. 
This is shown in Fig. |7| (a), where the probabilities of 
neighboring states, (n, sign{s)) = (n — 1,-1-) and (n, — ), 
oscillate at a frequency ^ 2y^g. In practice, however, 
each of those amplitudes is of course also weakly coupled 
to others due to finite tunneling. Note that the proba- 
bility of the (0, — ) state remains almost constant at all 
times since it has no direct partner. 

As zJb/g increases, the system dynamics becomes in- 
creasingly dominated by the molecular condensate com- 
ponent. In the extreme case of zJb/ g = 10, we reach a 
regime of clear periodic motion between the molecular 
and atomic superfluids. Figures [3(b) and (c) show that 
the Jaynes-Cummings character of Fig. [7| (a) disappears 
with increasing zJb/g^ and a high degree of coherence 
develops, accompanied by the formation of new order 
characterized by the simultaneous appearance of several 
molecular number states of same spin. 



B. Self-trapping transitions 

In this subsection we examine how the dissociation dy- 
namics is affected by the initial state of the molecular 
field, or stated another way, whether it undergoes a sud- 
den transition as a function of Ub/zJb- 

To answer this question we solved numerically the 
GW mean-field equations for several values of Ub/zJb 
for zJb/g — 20 and (5 = 0. Fig. |H1 shows the time evolu- 
tion of n;, (solid line), |4>p (dashed line), and |Ap (dot- 
dashed hne) from (a) Ub/zJb = 1.50 to (f) Ub/zJb — 
5.50. Since the initial superfluid fraction decreases as 
Ub/ zJb increases, these results illustrate situations where 
the dissociation dynamics changes from being "coher- 
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FIG. 6: (Color online) Time evolution of Uh (solid line), |<E>|'^ 
(dashed line), and |Ap (dot-dashed line) for (a) zJb/g = 0.1, 
(b) zJb/g — 1, and (c) zJb/g — 10, and for 5 — zJt, Ub/zJt — 
0.1. 




FIG. 7: Gray scale rendition of the time evolution of the 
probabilities |/n,sp of the states (n, sign(s)) for (a) zJt/g = 
0.1, (b) zJb = 1, and (c) zJb/g = 10, respectively. The darker 
shades of gray indicate higher probabilities. Same parameters 
are taken as in Fig. |S| 



ent", or condensate-dominated, to "incoherent," or non- 
condensate dominated. Interestingly, we find two types of 
transitions against small variations of Ub/zJi,: the first 
one is a sudden suppression of the amplitude of oscil- 
lations between Figs. |H1 (b), Ub/zJb — 1.86, and (c), 
Ub/zJi, = 1.87; and the second is a sudden increase of 
that amphtude between Figs.|Hl(d), Ub/zJb = 4.70, and 
(e), Ub/zJb = 4.71. Such a transition from a small ampli- 
tude oscillation (localized solution) to a large amplitude 
oscillation (delocalized solution), and vice versa, due to 
a nonlinearity of the system is known as a self-trapping 
transition [isl llq . Since these two transitions have very 
different coherence properties, we refer to them from now 
on as the coherent self-trapping transition for the former 
and the incoherent self-trapping transition for the latter, 
respectively. 

The Tavis-Cummings model of Sec. Ill provides a use- 
ful way to understand why two different self-trapping 



transitions are present against a variation of Ub/zJb- If 
we regard the third term in Eq. (|2DJ), corresponding to 
the molecule-molecule interaction, as a nonlinear detun- 
ing, the energy difference between number states of n 
condensate bosons and n — 1 condensate bosons is given 
by hfln — hSo + hUa{n/Ns) where we have neglected a 
term of order l/Ng- The dynamical properties are de- 
termined by the competition between that detuning and 
the photoassociation term in Eq. (^0)1 . When the detun- 
ing is dominant the dynamics is characterized by small 
amplitude oscillations, while large amplitude oscillations 
appear when the atom-molecule coupling is dominant. In 
the intermediate region where these two effects are com- 
parable, that is, for rj„ ~ a self-trapping transition 
appears. Replacing n/Ns in the nonlinear detuning en- 
ergy f2„ by the initial superfluid density of bosons per 
site, \^{t — 0;C/fc/zJf,)p which depends on the ratio of 
Ub/zJb as shown in Fig. |21 we obtain, for S = and 
C/fc/ =0, 

n„ = -zJb + UbMi = 0; Ub/zJb)\^ 

This frequency becomes equal to zero for two values of 
Ub/zJb- For large values of zJb/g, such as in Figs.|Hl self- 
trapping transitions occur around at these two resonance 
points. 

Self-trapping transitions can also be understood in 
a different picture from the semiclassical form H21|) of 
the generalized Tavis-Cummings Hamiltonian, which de- 
scribes the motion of a classical particle in the potential 
V{no) H22|) . At t = 0, the particle is at rest and is lo- 
cated at some point uq = \^{t — 0)p on the potential. 
Under the influence of the dissociation laser it then starts 
rolling down the potential y(no) toward smaller values 
of no, passes through a potential minimum, and moves 
uphill up to a point having the same potential energy as 
initially. At this point the particle turns back, result- 
ing in periodic oscillations in the (continuous) molecule 
number. What happens near the self-trapping condition 
is that an additional potential barrier appears in the po- 
tential. For delocalized (resp. localized) solutions, the 
particle can (resp. cannot) climb up the barrier, resulting 
in a sudden transition between large and small amplitude 
oscillations. The key factor in achieving this transition 
is the quartic form of the potential T^(no), which disap- 
pears in the absence of the effective two-body interaction 
Uo 40].) 

This classical interpretation relies on the assumption 
that the dynamics of the condensate fraction is uncou- 
pled to that of the non-condensate component. Since 
the GW variational state also involves a non-condensate 
component, the corresponding mean- field solutions are 
expected to deviate somewhat from those of Eq. 121|l . 
especially in the vicinity of the incoherent self-trapping 
critical point. Figure|51shows the critical curve of the self- 
trapping transitions in the Ub/zJb - zJb/g plane for (5 = 0. 
The dashed line corresponds to the semiclassical solution 
of the Tavis-Cummings dynamics, obtained by replac- 
ing J/ ^ j/Q = |<i>(t = 0)p in the potential V{no). The 
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FIG. 8: (Color online) Time evolution of n;, (solid line), [$|^ (dashed line), and lAI'^ (dot-dashed line) for zJb/g = 20, 5 = 0, 
and for (a) Ub/zJi, = 1.50 (b) Ub/zJb = 1.86, (c) Ub/zJb = 1.87, (d) Ub/zJb = 4.70, (e) Ub/zJb = 4.71, and (f) Ub/zJb = 5.50. 



stars are mean-field results obtained by seeking a point 
where a sudden suppression of amplitude happens in the 
first oscillation period Both these results have a C 
shape, with localized solutions present inside the curve, 
and share several common features: first, we observe that 
two transitions occur for a given zJb/ g. They correspond 
to coherent^ respectively incoherent transitions, for the 
lower and upper values of Uh/zJ\,. Second, the transi- 
tions disappear below some value of zJb/ g- 

In the vicinity of the incoherent transition, the mean- 
field dynamics of the delocalized solutions is very dif- 
ferent from its semiclassical counterpart. As shown in 
Fig. IHl^e) , a large incoherent oscillation occurs once the 
molecular condensate reaches a sufficiently small value, 
indicating that the delocalization in the motion of molec- 
ular condensate is a trigger for an incoherent large ampli- 
tude oscillation in the system. This region is also char- 
acterized by large fiuctuations 



n\Y) 



in the molecule number fiuctuations, a feature that is of 
course ignored in the semiclassical description. Fig. [TUI 
shows the time evolution of for the cases of Fig. |H1 (d) 
(dashed line), and (e) (sohd fine). 



C. Intersite spin-spin coupling 

Finally, we briefly discuss the effects of XXZ term in 
the Hamiltonian Hs- It is known that the trap frequen- 
cies of homonuclear molecules in a lattice site almost co- 
incide with the atomic trap frequencies 0, Il4j , so that in 
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FIG. 9; Critical curve of the self-trapping transitions. The 
star and dashed lines are for the GW mean-field and the semi- 
classical Tavis-Cummings solutions, respectively. The inside 
of the curve is the domain of localized solutions. 



general J/ ^ J& due to the higher mass of the molecules, 
TOfc = 2to/. However, if the inequality |?7/| 3> zJf that 
underlies all of the considerations in this paper is ful- 
filled, the effective tunneling of the fermions is given by 
Js = AJj/Uf, which can then easily lead to a system of 
\Js\ < Jb where the intersite spin-spin coupling does not 
affect the coherence properties of the dynamics. In the 
following we limit our considerations to this case 'S] . 

In the regime of cooperative dynamics of Fig. [S] the 
spin-spin coupling gives rise to a positive (resp. nega- 
tive) detuning effect for the case of attractive (resp. re- 
pulsive) fermion-fermion interaction, since a superfiuid 
state of Fermi pairs (| A| 7^ 0) lowers, resp. raises the en- 
ergy for Uf < 0, resp. Uf > 0. This additional detuning 
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FIG. 10: Molecule number fluctuations per site, ai,, across 
the incoherent self-trapping transition point for the cases of 
Fig. |H| (d) (dashed line) and (e) (solid line) . 

can however be compensated by changing the photodis- 
sociation detuning energy 6. We have checked that finite 
vahies of Js preserve the appearance of the self-trapping 
transitions. Finally, in the regime of incoherent dynamics 
the effects of the XXZ term become weaker with decreas- 
ing zJb/g. We thus conclude that the discussion in the 
preceding subsections holds even in the presence of the 
intersite spin-spin coupling. 



two-component fermionic atoms in an optical lattice. 
We showed that the dissociation dynamics in the strong 
fermion-fermion interaction limit, \Uf\ ^ zJ/, is gov- 
erned by a spin-boson lattice Hamiltonian dynamics. Us- 
ing a Gutzwiller mean-field analysis, we examined the 
crossover of the dissociation dynamics from an indepen- 
dent single-site regime to a cooperative regime as the 
molecular tunneling is increased. We also showed the 
appearance of two types of transitions, that is coher- 
ent and incoherent self-trapping transitions, as a small 
variation of Ub/zJb- The mean-field approximation only 
accounts for those tunnel processes proportional to the 
expectation values of bosonic and fermionic pair opera- 
tors. Thus, when starting from the MI molecular phase 
the dynamics of the atoms and molecules inside the in- 
dividual lattice sites always remains independent of the 
other sites. In this situation, genuine quantum effects 
are significant, and are closely related to the stability of 
the MI phase under boson-fermion conversion 19] and to 
the formation of SF correlations across the MI-SF critical 
point This will be discussed in a future publication. 
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